arXiv: 1507.00425vl [astro-ph.GA] 2Jul2015 


Astronomy & Astrophysics manuscript no. CaTcal'vl4 © ESO 2015 

July 3, 2015 


The Calcium Triplet metallicity calibration for galactic bulge stars. * 

S. Vasquez 12 , M. Zoccali 1,2 , V. Hill 3 , O. A. Gonzalez 4 ,1. Saviane 4 , M. Rejkuba 5 , and G. Battaglia 6 ' 7 

1 Instituto de Astroffsica, Pontificia Universidad Catolica de Chile, Av. Vicuna Mackenna 4860, 782-0436 Macul, Santiago, Chile 
e-mail: svasquez@astro. puc. cl 

2 Millennium Institute of Astrophysics, Av. Vicuna Mackenna 4860, 782-0436 Macul, Santiago, Chile 

3 Laboratoire Lagrange (UMR7293), Universite de Nice Sophia Antipolis, CNRS, Observatoire de la Cote d’Azur, CS34229, 06304, 

Nice Cedex 04, France 

4 European Southern Observatory, Av. Alonso de Cordova 3107, Casilla 19, 19001, Santiago, Chile 

5 European Southern Observatory, Karl-Schwarzschild Strasse 2, D-85748 Garching, Germany 

6 Instituto de Astroffsica de Canarias, calle via Lactea s/n, 38205 San Cristobal de La Laguna, Tenerife, Spain 

7 Universidad de La Laguna, Dpto. Astroffsica, 38206 La Laguna, Tenerife, Spain 

Preprint online version: July 3, 2015 


ABSTRACT 

Aims. We present a new calibration of the Calcium II Triplet equivalent widths versus [Fe/H], constructed upon K giant stars in the 
Galactic bulge. This calibration will be used to derive iron abundances for the targets of the GIBS survey, and in general it is especially 
suited for solar and supersolar metallicity giants, typical of external massive galaxies. 

Methods. About 150 bulge K giants were observed with the GIRAFFE spectrograph at VLT, both at resolution R~20,000 and at 
R~6,000. In the first case, the spectra allowed us to perform direct determination of Fe abundances from several unblended Fe lines, 
deriving what we call here high resolution [Fe/H] measurements. The low resolution spectra allowed us to measure equivalent widths 
of the two strongest lines of the near infrared Calcium II triplet at 8542 and 8662 A. 

Results. By comparing the two measurements we derived a relation between Calcium equivalent widths and [Fe/H] that is linear over 
the metallicity range probed here, -1 <[Fe/H]< +0.7. By adding a small second order correction, based on literature globular cluster 
data, we derived the unique calibration equation [Fe/H] Ca T = -3.150 + 0.432 W + 0.006IF' 2 , with a rms dispersion of 0.197 dex, valid 
across the whole metallicity range -2.3 <[Fe/H]< +0.7. 

Key words. Stars: abundances - Galaxy: bulge - Techniques: spectroscopic 


1. Introduction 

The Calcium II Triplet (CaT) at ~8500 A is one of the most 
widely used metallicity index, as well as an excellent feature 
to measure radial velocity at low spectral resolution. The three 
lines at 48498, 48542, 48662 A are so strong that they can be 
measured easily at low resolution and at relatively low signal 
to noise. In addition, their location in the near-infrared part of 
the spectrum is ideal to observe the brightest stars of any not 
too young stellar population, namely cool giants. CaT spectra 
of cool giants can be obtained with reasonable exposure time 
both for external galaxies in the local group, too far away to be 
observed at high spectral resolution, and for Milky Way stars in 
high extinction regions, such as the Galactic bulge. 

Obviously the popularity of the CaT spectral feature re- 
sides in how accura t ely it can be used to measure metallicities. 
lArmandroff & ZinTib 1 988/ first demonstrated that the equivalent 
widths (EWs) of CaT lines, in the integrated spectra of glob¬ 
ular clusters (GCs), strongly corre lated with the cluster m etal- 
licity [Fe/H]. A few years la ter, lOlszewski et al.l (Il99ll) and 
Armandroff & Da Costal (1199 ll) analyzed the behaviour of CaT 
lines in individual cluster stars. They noticed that the EWs of 
CaT lines show a dependence not only on metallicity, but also 
on absolute magnitude. They therefore introduced the use of “re¬ 
duced equivalenth widths” (W'), which corresponds to the sum 


* Based on observations taken with ESO telescopes at the La Silla 
Paranal Observatory under programme ID 385.B-0735(B). 


of some combination of the individual EWs, weighed by the dif¬ 
ference between the star V magnitude and the magnitude of the 
Horizontal Branch in the same cluster (V - Vhb)- 

Several empirical relations between the reduced EWs of CaT 
lines and the [Fe/H] abundance are present in the literature. Most 
of them used star clusters, for which [Fe/H] abundance could be 
derived in several ways, and not necessarily for the same stars for 
which CaT was measured. Among those, a very comprehensive 
one is the work of lRutledge et al.l (119971) who derived CaT metal¬ 
licities for 52 Galactic GCs in a homogeneous scale, covering 
the range -2 <[F e/H]< -0.7. By comparing their scale with the 
classical lZinn & Westl d 1984 ) metallicity scale, they found a non¬ 
linear relation between the two. Conversely,_a linear relation was 
found between the CaT metallici ties by Rutledge et alj ( 19971) 
and the metallicities derived bv lCarretta & Grattonl (1991 ). The 
latter are based on Fe lines, measured on high resolution spectra. 

Traditionally, the CaT metallicity calibration has been con¬ 
structed based on RGB stars in globular (Cole et al. 2004: 


Warren & Co 

e 2009; Saviane et al l 2012; 

(Carrera et al 

2007, 201^"Open clusters 


to extend the metallicity range up to [Fe/H]—1-0.5, at the same 
time probing a younger age regime (13 Gyr<age<0.25 Gyr). 
The CaT EWs seem to be only weakly dependent on the age of 
the star, although only a few star clusters constrain the relation 
at high metallicity, and anyway old star clusters at supersolar 
metallicities are not available to set a robust constrain on the age 
dependance. 
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[Battaglia et al.l (2008) compared a CaT calibration based on 
GCs with one based on red giants in dwarf spheroidal galax¬ 
ies. The latter are complex stellar populations with a range of 
ages and, most importantly, a range of [Ca/Fe] abundance ra¬ 
tio. They conclude that the two relations are fairly consistent 
within the errors, although small differences exist. In particu¬ 
lar, a CaT-to-[Fe/H] calibration based on GCs would overesti¬ 
mate the [Fe/H] of low metallicity stars ([Fe/H]< -2.2) by -0.1 
dex, while underestimating the [Fe/H] of (relatively) high metal¬ 
licity stars ([Fe/H]> -1.2) by -0.1-0.2 dex. It is important to 
emphasize that their analysis is limited to the metallicity range 
-2.5 <[Fe/H]< -0.5. 

More recently [Starkenburg et al.l (i2010t) presented a syn¬ 
thetic spectral analysis of the CaT method to derive iron abun¬ 
dance, with particular emphasis at the low metallicity regime. 
They derive a new calibration, valid in the range -4 < [Fe/H] < 
-0.5, that deviates from the commonly used linear relation for 
[Fe/H]< -2. 

To summarize, CaT to [Fe/H] calibrations at solar and super¬ 
solar metallicities are very sparse, and based only of young open 
clusters. Deviations from linearity, in the relation between CaT 
EWs and [Fe/H] in the high metallicity regime would not be sur¬ 
prising, given that the Calcium lines are highly saturated already 
well below solar metallicity. Departures from local thermody¬ 
namical equilibrium are also known to affect the core of the lines 
(iMashonkina et al.ll2007l) because of saturation. In this paper we 
probe the high metallicity regime of this relation based on old 
red giants in the Galactic bulge. 

This study is the result of a pilot project we carried out to 
investigate the possibility of deriving the metallicity distribution 
function (MDF) of the Galactic bulge, based on an iron abun¬ 
dance scale, by measuring CaT lines in low/intermediate res¬ 
olution GIRAFFE spectra. To this end, we acquired both high 
(R~20,000) and low (R~6,000) resolutions GIRAFFE spectra 
for a large sample of K giants in two fields in the Galactic bulge. 
When preliminary results were very encouraging, we were 
granted an ESO Large Programme (187.B-0909; PI: Zoccali) to 
observe 200-400 red clump (RC) giants in each of 25 fields in the 
Galactic bulge, with GIRAFFE at low resolution. Those were 
then combined with other archive data obtained with the same 
instrument and setup, also targetting bulge RC stars. The result¬ 
ing dataset, including spectra for 6392 stars in 31 fields, is what 
we called the Giraffe Inner Bulge Survey (GIBS), described in 
JZoccali et al.ll2014l hereafter Paper I). Paper I also lists the ra¬ 
dial velocity of individual stars in the sample, and presents radial 
velocity and velocity dispersion maps of the inner bulge. 

We focus here on the CaT to [Fe/H] calibration, that will be 
used in a forthcoming paper to derive the MDF in 25 fields in 
the Galactic bulge. 


2. Observation and data reduction 

The target stars we used to construct the CaT calibration are all 
located in Baade’s Window, and they consist of two samples: 
80 stars in the bulge RC, and 68 stars on the RGB, -0.7 mag 
brigther than the RC. The medium resolution spectra analysed 
here were taken with the FLAMES-GIRAFFE multifibre spec¬ 
trograph at the VLT, in Medusa mode, within the ESO program 
ID 385.B-0735(B). The spectrograph allows us to obtain spectra 
for -130 objects, including sky positions, over a 25 arcmin di¬ 
ameter field of view. We selected the LR08 setup, which covers 
the wavelength range between 8206 and 9400 A, with a spectral 
resolution of R « 6500. The exposure time needed to reach a 



v-i 

Fig. 1 . The target selection is shown over the V and I band colour 
magnitude diagram using color dots. Red and blues correspond 
to the RC and RGB stars from Baade’s window (b = -4). 


mean signal to noise of - 50/pix were 900 and 1800 seconds for 
RGB and RC stars, respectively. 

High resolution spectra for the targets in Baade’s Window 
were obtained from Zoccali et al. (2008) and Hill et al. (2011) 
who used the same instrument in the setups HR13 and HR14, at 
resolution R - 22,500 and R - 28,800, respectively. Figure Q] 
shows the corresponding color magnitude diagrams (CMDs) for 
the target selection. 

The data reduction was done using the instrument pipeline 
provided by ESO (version 2.8), which performs bias, flat field 
corrections, wavelength calibration, and 1-D extraction based on 
the daytime calibration images. Since the pipeline does not in¬ 
clude the sky subtraction routine, we performed this task with 
IRAF as follows. As a first step, a master sky spectrum was 
created by combining all the sky spectra (~ 20 per setup) us¬ 
ing the median and a sigma clipping algorithm. Figure [2] (up¬ 
per panel) shows the raw spectrum of a typical star, around 
the CaT region, together with the master sky spectrum for the 
same exposure. It is clear that the three CaT lines do not suffer 
from strong contamination from sky emission lines compared 
with CaT line strength. The master sky spectrum was then sub¬ 
tracted from each individual science spectrum using the IRAF 
task skytweak. The latter allows shifting and scaling of the mas¬ 
ter sky spectrum until a good subtraction is achieved. In order 
to find the best shift and scale factor we used the strong emis¬ 
sion lines located in the wavelength range between 8261 to 8481 
A and 8739 to 9071 A. 

The spectra were normalized with the IRAF task 
continuum, using the wavelength range between 8250- 
8461 A, 8568-8650 A, and 8700-9273 A, in order to avoid the 
CaT lines. Finally, radial velocities were measured by means 
of cross-correlation (fxcor task) against a template synthetic 
spectrum for a typical metal poor K giant star (T e g = 4750, 
log g = 2.5 and [Fe/H] = -1.3). Such low metallicity was 
chosen in order to ensure that the CaT lines in the template were 
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Fig. 2. The upper panel shows the normalized raw spectrum for 
a typical Milky Way bulge star, together with the corresponding 
sky spectrum. The lower panel shows processed spectra for 4 
bulge stars spanning the typical metallicity range of the Galactic 
bulge. The increasing strength of the CaT lines with metallicity 
is evident, together with the increasing number of small atomic 
lines contaminating the continuum. 


Fig. 3. The upper panel shows for comparison the more metal 
poor and metal rich synthetic spectra from our library, in the 
region around the CaT lines. In both cases, the black line shows 
the spectra including the lines of all elements, while the red lines 
shows only the three CaT lines. The lower panels show the dif¬ 
ference between the real EW, as measured in the red spectra, and 
the measured EW, as measured in the black spectra, for each of 
the three CaT lines, as a function of EW and for three values of 
S/N. 


relatively deep but free from contaminations by other lines on 
their wings. 

Example of the final spectra, sky subtracted, radial veloc¬ 
ity corrected and continuum normalized are shown in the lower 
panel of Fig. [2] for stars of different metallicities. An arbitrary 
shift has been applied along the y-axis, to avoid overlaps. 


3. The Calcium Triplet calibration 

The spectral region around the CaT feature (8350-9000 A) 
includes several lines which might dominate the spectrum, 
depending on the spectral type of the star under analysis. 
ICenarro et al.l fePOlil showed that stars with spectral type be¬ 
tween F5 (T eff ~ 6500 K) and M2 (7’ e ir ~ 3700 K) are the best 
targets to measure the CaT lines, due to their only marginal con¬ 
tamination by Fel, Mgl and Til lines. CaT lines in hotter stars are 
blended with Paschen lines, while in cooler stars the continuum 
is plagued by TiO molecular bands. Being of K spectral type, 
our targets are well inside the recommended temperature range. 

Additional constraints have been proposed in terms of the 
brightness range where a si ngle gravity c orrection can be used 
(see section 3.2 for details'). [Carrera et alJ ( 2007 1 defined a min- 
i mum abso lute magnitude My < 1.25 or Mj < 0.0, while 
iDa Costa et al.l d2009t) recommends a lower brightness limit of 
~ 0.2 mag below the RC (or the horizontal branch). Both our 
RGB and RC targets are inside these constraints, with the RC 
stars being within the intrinsic dispersion - mainly due to metal¬ 
licity and distance spread - of the RC in Baade’s Window. 


3.1. The equivalent width measurements 

While both our RC and RGB target stars are inside the recom¬ 
mended temperature and gravity range for CaT measurements, 
the metallicity distribution of the Galactic bulge reaches (ex¬ 
tends) to significantly higher metallicity than typical star clus¬ 
ters for which this calibration has been defined and used so far. 
The impact of the metal content over the CaT EWs measure¬ 
ments is mostly driven by the continuum depletion due to the 
increasing contamination of the CaT lines by other atomic lines. 
In order to explore the impact of the metal content on the CaT 
EWs measurement, we analysed a grid of synthetic spectra gen¬ 
erated using turbo spec trum. The latter is a code of line synthe¬ 
sis described in I Alvarez & Plezl (1998) and improved along the 
years by B. Plez. In this case it was used together with the grid 
of OS MARCS spherical model atmosphere^] (Gustafsson et al. 
2008). The synthetic spectra have the same resolution as the 
GIRAFFE spectra, and a set of stellar parameters representative 
of a K giant star (7' c n = 4750 and log g = 2.5). Only the metal¬ 
licity was allowed to vary across the range from [Fe/H]=-1.5 to 
the extreme case of [Fe/H]=+1.0 detfl keeping [Ca/H]=[Fe/H] 
in the whole grid. 

Synthetic spectra for stars at the two extreme ends of the 
metallicity grid are plotted in the upper panel of Fig0 The black 
line is the synthetic spectrum containing list of all the element in 
the wavelength range between 8450 and 8700 A. The red line, 

1 models available at http://marcs.astro.uu.se/ 

2 The metallicity distribution of the Milky Way bulge spans the range 
between -1.5 up to +0.7 dex 
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Table 1. 


Feature (A) Line bandpass (A) 

8498.02 8491.02-8505.02 

8542.09 8530.09 - 8554.09 

8662.14 8653.14-8671.14 


in contrast, shows only the three CaT lines. This exercise clearly 
shows the depletion of the continuum (pseudo-continuum), at 
high metallicity, due to the presence of many unresolved lines. 
Because we deal with normalized spectra (without flux calibra¬ 
tion) we need to include the effect of continuum depletion in the 
empirical calibration itself. 

The measurements of the CaT EWs can be done by fit¬ 
ting an empirical function to each line profile and calculate 
the corresponding EW from direct integration of the func¬ 
tion. Different functions have been explored over the years. 
[Armandroff & Da Costal (119911) used _a Gaussian profile, when 
fitting stars in globular clusters. [Rutledge et al.l (1997), instead, 
used a Moffat function. Both functions give good results only 
for stars of relatively low metallicity. In fact, ICole et all d2004l) 
showed that in the metal rich regime, the strong wings of the 
CaT lines can be reproduced by the sum of a Gaussian com¬ 
ponent, fitting the line core, and a_ Lorentzian function, fitting 
the wings. lErdelvi-Mendes & Barbuvl ( 1991 ) and lBattaglia et al.l 
( 2008 ) showed that the EWs of CaT lines is driven mostly by 
the wings, because of the strong pressure broadening. In other 
words, the strength of the line is mostly driven by the electronic 
pressure (i.e., metallicity), rather than by the usual combination 
of temperature and elemental abundance, shaping weak lines. 

In order to test the accuracy of the Gaussian+Lorentzian fit 
across a range of metallicity and signal to noise, we performed 
some measurements on the grid of synthetic spectra described 
above. As a first step, EWs were measured in the synthetic spec¬ 
tra containing only Calcium lines by means of numerical integra¬ 
tion. These were assumed to be the real EWs. As a second step, 
the grid of synthetic spectra was convolved to the resolution and 
resampled in wavelength to match real data and degraded to dif¬ 
ferent S/N, from 10 to 100 in steps of 10. A Poisson noise model 
was used for the last operation. The resulting spectrum was then 
re-normalized to a pseudo-continuum of 1, and the combina¬ 
tion of a Gaussian and a Lorentzian function was fitted to the 
CaT l ines. The wavelength ranges used by lArmandroff & Zinnl 
d!988l) for the fit were slightly modified, as listed in Table Q] in 
order to improve the fit of the wings. The lower panel of Fig. [3] 
shows the difference between the measured and the real EWs, 
for each of the CaT lines, as a function of EW (metallicity) and 
S/N. It is clear that the EWs of the two strong lines would be un¬ 
derestimated in the data, and more so at high metallicities. This 
systematic reaches up to ~1 A around [Fe/H]=+0.5. However, 
the trend is very smooth and it is not affected by S/N, down to 
a S/N-10, much lower than our poorest spectrum. As a con¬ 
sequence, a correction for this systematic will be implicitly in¬ 
cluded in the [Fe/H]-vs-CaT empirical calibration equation de¬ 
rived in Section 3.4. For the weakest line, the systematic differ¬ 
ence between real and measured EW is much smaller than for 
the strong lines. However, a step is found around EWj = 1.5 A, 
due to the appearance of some blends in the wings of this line. 
Additionally, the weakest line is the most affected by unperfect 
sky subtraction, due to the presence of two sky lines close to its 
wings. For this reason, we discarded this line in the following 
analysis. 



Fig. 4. Fitting examples for three bulge stars covering the metal¬ 
licity range of the Milky Way bulge. The red line correspond 
to the best fitting function found for each line, while as dashed 
lines are plotted the individual contribution from Gaussian and 
Lorentzian functions. 



Fig. 5. The strength of the two strongest CaT lines versus the 
sum of the two lines, X EW. The corresponding fractional con¬ 
tribution with respect to X EW is shown in each panel. 


Examples of line fitting are shown in Fig. |4] for three ob¬ 
served spectra covering the metallicity range of the Milky Way 
bulge. For each star, the big black dots show the observed spec¬ 
trum, the red line is the final fit, while the thin dotted and dashed 
lines show the individual contribution of the Gaussian and the 
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SEW = EW„, 


EW 8682 (A) 


gravity of the star. Both lArmandroff & Da Costal (1991) and 
[Olszewski et al.l f 1 99 lb showed that the combined effect of ef¬ 
fective temperature and surface gravity on the EWs of CaT 
lines can be parametrized as a function of luminosity alone. 
The CaT EWs corrected by this effect (i.e., the reduced EWs) 
depend only upon metallicity. Specifically, they found that, in 
each given cluster, red giant stars follow a straight line in the 
luminosity-IEW plane. The slope of the line is identical for dif¬ 
ferent clusters, while its zero point depends on the cluster metal¬ 
licity. As a proxy for lumin osity, some authors use the abso¬ 
lute magnitude of the star (iCarrera et al.ll2007l 1201 3b . or, most 
commonly, the difference between the star apparent magnitude 
and the magnitude of the cluster horizontal branch, o r red clump 
( ArmandrcS^&Da_Costa| j~991j Rutledge_et al] 1997 : Cole_etal 


2004 : Battaglia et al. 2008 : IWarren & Colei 2009 Is aviane et al 


Fig. 6. Comparison between our measurements and the literature 
values of the sum of the two stronge s t lines for different stel¬ 
lar clu ster taken from IWarren & Colei (12009b and ISaviane et al.l 

d20l2h . 


Lorentzian profiles, respectively. Additionally, the line-ratio be¬ 
tween the two strongest CaT lines is remarkably constant across 
the metallicity range spanned by the bulge data. Figure [5] shows 
the contribution of the two strongest individual CaT lines to the 
sum of the two. The data can be fitted by a single line, with the 
slope indicated in the label. This line ratio can be used, for in¬ 
stance, to predict the EW of one of the lines, when the other is 
not available for a direct measurement. It can also be used to 
test the consistency of the measurements, in order to reject lines 
with anomalous line ratio. This is particularly useful when deal¬ 
ing with large surveys. 

As will be explained in the following sections, we will use 
cluster stars to constrain the final metallicity CaT calibration. 
Before doing it, we need to verify that our EWs measurements 
are in the same scale as those used in previous works. To this 
end, we obtained the spectra of the two most metal rich clusters 
(NGC 6791 and NGC 6819) in the sample bv IWarren & Colei 
(12009b and re-measure EWs. The spectra kindly provided to 
us, by private communication, were already wave-calibrated 
and sky-subtracted, therefore we only applied the same contin¬ 
uum normalization procedure and EW measurement used for 
our data. Additionally, in order to to probe a wider range of 
metallicity, we also measured a sample of globular cluster stars, 
drawn from lSaviane et al.l d2012b . spanning the metallicity range 
-2.3 <[Fe/H]< +0.07. The selected stars belong to the globu¬ 
lar clusters M 15, NGC 6397, M4, NGC 6553, and NGC 6528. 
Figure [6] shows the excellent agreement between our measure¬ 
ments and the literature one, for the sum of the EWs of the two 
strongest CaT lines, hereafter SEW (= EW ^542 + EW^ei)- 

3.2. The reduced equivalent width 

The strength of an atomic line, in general, is driven by the el¬ 
emental abundance, the effective temperature, and the surface 


2012 : Mauro et al. 2014). We adopted the second option, in or¬ 
der to minimize the uncertainty in the distance and reddening, 
both variable across the Galactic bulge. 

The K s magnitude, from the VVV survey (Minniti et all 
1201 ob . was used for our bulge stars. In order to determine the 
slope in the luminosity-XEW plane we used the star cluster data 
from I Warren & Cold (2009). For each cluster star we plot FEW 
as a function of (A- Arc) and fit a linear relation, deriving a rep¬ 
resentative slope by making the average of the single measure¬ 
ments. The left panel of Fig. [7] shows the straight lines defined 
by different clusters, each in a different color. They are consistent 
with a constant slope of J3k s = 0.384 + 0.019 A ma gT 1 , in ex¬ 
cellent agreement wi th those derived by Mauro et all (12014b us¬ 
ing the clusters from lSaviane et all (l2012h combined with VVV 
photometry. This slope allows us to define the reduced EW pa¬ 
rameter 

W' = SEW + Pk s (K - Arc) 

which is dependent only upon the metallicity of the star. 

In order to determine the corresponding (A - Arc) for our 
bulge stars we measure the mean RC magnitude from the com¬ 
plete VVV de-reddened K s photometric catalogue correspond¬ 
ing to the observed field, by using the reddening map presented 
in lGonzalez et alJ ( 201 l c). The right panel of Fig. [7] shows the 
location of our RC and RGB targets, all in Baade’s Window, 
in the same plane. The color coding traces metallicity (as mea¬ 
sured from the HR spectra), and it is the same used for the star 
clusters. Contrary to what happens in star clusters, bulge stars in 
each metallicity bin do not seem to define parallel straight lines. 
We interpret this as the result of the probable disk contamina¬ 
tion, the errorbar on the metallicity of each star, and the intrinsic 
width of the metallicity bin. Also, while for star clusters the x- 
axis is a clean proxy for luminosity, bulge stars have a significant 
distance spread, and the red clump only reflects their mean dis¬ 
tance. 


3.3. Reference metallicities: Galactic bulge and globular 
clusters 

In order to derive an equation converting between W and 
the [Fe/H] abundances, a set of independent determinations of 
[Fe/H] from high resolution spectra are needed. In what fol¬ 
lows we describe how those were obtained, both for bulge field 
giants_and for a su bset of the globular clusters presented in 
IWarren & Cold d2009i) . The latter have been added to our sam¬ 
ple, after showing that they were fully compatible, in order to 
extend the metallicity range on the metal poor side. 
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Fig. 7. The CaT equivalent widths of the star cluster RGBs (left panel) and Galactic bulge stars (right panel) are compared with their 
corresponding magnitude above the mean red clump magnitude. The colour code is used to mark clusters of different metallicity, 
together with their best-fitting lines, assuming the best common slope of 0.384 Afhag -1 . Thin dotted lines in each panel show this 
slope, to guide the eye. 


3.3.1. Homogeneous high resolution metallicities 
determination for Galactic bulge stars 

[Fe/H] abundances for both RC and RGB sta rs i n Baade’s 
Windo w, were discussed in IZoccali et alJ ([2008) and iHill et alJ 
(1201 ll) . They were based on GIRAFFE spectra, through setups 
HR13 and HR14, at resolution R~22,500 and ~28,800, respec¬ 
tively. After the pre-reduction and coaddition of the spectra ex¬ 
plained in Section 2, iron abundances were measured by means 
of individual, unblended iron lines. Conceptually, we used the 
same method explained in our previous works, cited above. 
Namely, first guess effective temperature and gravity were de¬ 
rived from photometry, and used to construct an approximate 
model atmosphere. The latter was used to derive [Fe/H] abun¬ 
dances, independently, from several Fel and Fell lines. The 
stellar surface parameters (T e ff, log g and the microturbulent 
velocity x) were then refined by imposing consistency among 
the abu ndances derived from different lines. As discussed in 
iHill et all (1201 ll) . the actual implementation of this method may 
yield slightly different final abundances -especially for metal 
rich, cold giants- depending on which of the stellar parameters 
is fixed first. In other words, there can be degeneracies in the pa¬ 
rameter space, due to correlations among the stellar parameters. 

In order to minimize this effect, for the present work we 
applied an_automated procedure, based on the use of GALA 
dMucciarelli et al. i 12013 1). This code is designed to impose the 
excitation and ionization equilibria on Fel lines, thus deriving 
T e ff and log g, respectively, and to converge on a value for the 
microtubulent velocity that yields the same [Fe/ Hl ab undance 
for lines of different EWs. The ATLAS9 ( Castelli & Kurucj 
l2004l> model atmospher e were used instead of the MARCS ones 
iGustafsson et al.ll2008b because GALA is optimized for the for- 
mer. T he re sults are fully compatib le with those of lZoccali et al.l 
(120081) and iGonzalez et al.l (1201 lch . A small offset -within the 


errorbars- is present between the resulting metallicities and those 
presented bv IHill et al.l (1201 ll) . The use of GALA, however, does 
reduce at least the statistical errors, as demonstrated by the ex¬ 
ceptionally small spread of abundances in the [Mg/Fe] versus 
[Fe/H] plot, presented in IZoccali et al.l d2015l ) and Gonzalez et 
al. (2015, in preparation). 


3.3.2. Globular clusters 


In order to constrain the CaT versus [Fe/H] calibration in the 
metal poor regime, where bulge stars are sparse, we considered 
adding to our sample some globular clusters. A natural choice 
was a subset of the [Warren & Colei ( 2009 1 cluster sample. This 
is for two main reasons: i) we had access to their spectra, and 
indeed we have shown in Fig. [6] that the EWs measured with 
our method is fully consistent with theirs; and ii) K s photome¬ 
try is available for all those stars. In their work. Warren & Colei 
( 2009 ) derived a calibration equation between the CaT reduced 
EWs, as der ived using the three lines, and [Fe/H] measured in 
the scale of ICarretta & Gratton (i 19971) . The latter scale, how¬ 
ever, has been revised in Carretta et al.l (12009b (hereafter C09), 
to what has now become a new standard (e.g. Harrisl [l996. 2010 
edition). We therefore updated the [Fe/ Hl ab undances to the C09 
scale using the compilation in ICarrera et al.l (120131) . For globu¬ 
lar clusters without new [Fe/H] measurements, we used instead 
the relation provided in C09, to transform the metallicity from 
the old to the new scale. It is worth emphasizing here that both 
globular clusters and bulge stars have the same calcium over 
iron ratio ([Ca/Eel—1-0.4) at [Fe/Hl~ -1 where they overlap 
(IGonzalez et al Jl201 1 at ICarrera et al.ll2013h . 

The calibration by iWarren & Colei ( 2009 ) includes several 
open clusters, constraining the metal rich regime. We excluded 
those from our analysis, because we sample the metal rich side 
with a large number of bulge stars, and our goal is to optimize 
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the calibration for bulge-like giants. In section 4 we will compare 
our calibration with open and globular clustr data. 

3.4. The metallicity calibration 

The calibration equation between W' and [Fe/H]HR was derived 
by means of an error-weighed least square fitting of the data 
points relative to individual stars, both in the bulge and in glob¬ 
ular clusters. This is different from what is usually done for star 
clusters, where the mean values of (W') and ([Fe/H]) for each 
cluster are used. Several fitting functions were tested before con¬ 
cluding that a linear relation was the best choice, because it min¬ 
imizes the artificial saturation of the calibration equation at high 
metallicity (i. e. a constant metallicity at higher W'). The lat¬ 
ter, indeed, may create an unphysically sharp cut-off in the iron 
distribution function. 

As a first step, a linear fitting was performed between W' and 
[Fe/H]HR for RGB and RC stars, independently. The results are 
shown in the left panels of Fig. [8] together with the rms (-0.23) 
which is comparable to the uncertainties of the [Fe/H]HR mea¬ 
surements. The fits are consistent with each other, and can there¬ 
fore be combined in the single linear relation shown in the right 
panel of the same figure, namely: 

[Fe/H] RGB+ Rc = (-3.504 ± 0.147) + (0.524 ± 0.023)W' (1) 

whith a rms of 0.233 dex. 

As a second step, 7 globular clusters were added to the plot, 
covering the metallicity range -2.3 <[Fe/H]< -0.7. The fit of 
a linear relation between the mean W' and mean [Fe/H], for the 
clusters only, yields: 

[Fe/H] GC = (-3.175 + 0.094) + (0.462 ± 0.025)W' (2) 

which, though slightly different from the equation found for the 
bulge (Eq.[]}, shows a smooth transition between the two when 
they overlap. This suggests that a single calibration can be used, 
valid across the whole range -2.3 <[Fe/H]< +0.7. Nevertheless, 
as globular cluster and Galactic bulge stars do not have the same 
data sampling (i.e. for globular clusters ( W ') and <[Fe/H | HR > 
are used instead of the measurements for individual stars) we 
cannot directly combine the data sets in order to fit a single 
function, because the resulting fit would be dominated by the 
metal poor data, having smaller errors in [Fe/H]. Because of 
this, we decided to create two evenly sampled fiducial data sets 
taken from the equations <jT]) and ([2]). in the W' range between 
0 < W' < 5.4 and 5.2 < W' < 10 respectively. These data sets 
were then combined into a single one and fitted by a second or¬ 
der function in order to take into account the small change of 
slope between the two linear equations (Eq.Q]and[2]). 

The resulting fit gives: 

[Fe/H] Ca T = -3.150 + 0.432W' + 0.006W' 2 (3) 

with a rms dispersion around the fit of 0.197 dex. No errors are 
quoted for the coefficients here, because, due to the procedure 
adopted to combine globular clusters and bulge stars, this is more 
a fiducial relation, characterizing the locus of the data in the W' 
versus [Fe/H] plane, rather than a formal fit. 

When applying the calibration, errors on the derived [Fe/H] 
values will be calculated as the quadratic sum of the errors on 
the EW, propagated through equation© and the rms of the same 
equation. Figure|9]summarizes our results, showing the final CaT 
vs [Fe/H] calibration, as a green line, together with the bulge 
stars (black dots) and the globular clusters (orange dots). For the 
latter, individual W' -measured by IWarren & Cole! ( 2009 )- are 



w 1 (A) 


Fig. 9. The CaT calibration for the combination between bulge 
stars (black) with globular clusters from [Warren & Colei (2009i) 
(orange). The upper left panel shows the corresponding calibra¬ 
tion metallicity residuals for bulge stars. 


plotted for individual stars, all at the same cluster <[Fe/H) IIK . 
The upper-left corner shows the histogram of the residuals for 
bulge stars, i.e., the differences between [Fc/H]hr and [Fe/H] Ga T- 
The distribution is symmetric around zero and the dispersion is 
consistent with the rms quoted above. 


4. Comparison with previous calibrations 

Different CaT vs [Fe/H] calibration have been proposed in the 
past, differing among each others mainly because of the differ¬ 
ent methods to measure EWs, and the different [Fe/H ] scales. A 
recent calibration is presented in lSaviane et al.l (12012)) who used 
14 globular clusters to define a cubic equation covering the range 
-2.33 <[Fe/H]< 0.07, in the C09 scale. Since the EWs mea¬ 
sured by us are fully compatible with the ISaviane et alJ (120121 ) 
(c.f.. Fig.©, the two calibrations can also be directly compared. 

Figure [TO] is color coded as Fig. [9] with bulge giants as 
small black dots and globular and open cluster mean values 
as orange and red dots, respectively ( Warren & (ole 120091) . 
Our calibration is shown as a g reen line. Also shown are the 
globular clusters in (ISaviane et all 12012) blue stars) and their 
calibration equation, as a blue dashed line. The two calibra¬ 
tions, as well as the original one from IWarren & Colei d2009i) . 
are consistent with each other for [Fe/H]< -0.8 dex. The 
higher order fit in ISaviane et al.l ( 2012 ) is required only by the 
three most metal rich globular clusters in their study, namely: 
NGC 5927 ([Fe/H]=-0.29), NGC 6553 ([Fe/H]=-0.16) and 
NGC 6528 ([Fe/H]=+0.07). Interestingly, the same behaviour 
is seen in C09, where the updated lRutledge et al] (119971) calibra¬ 
tion changes from a linear to a cubic fit when NGC 6553 is added 
to the original sample (see Fig. A.2 from lCarretta et all ( 2009 )). 
Moreover, as noted by Saviane et al., NGC 6528 has W' lower 
than NGC 6553, though being more metal rich. 
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Fig. 8. High resolution metallicities are compared with reduced equivalent widths of RGB (blue), RC (red) and the combination of 
the two (black) bulge star samples. The best fitting CaT calibration for each sample is also shown in each panel. 


Although we are not able to explain the difference be¬ 
tween the three most metal rich globular clusters and bulge 
stars, we note that the cubic relation of ISaviane et al.l d2012l) is 
clearly not appropriate for bulge stars, which are the main fo¬ 
cus of our interest. Also^ shown in Fig. [TO] are the open clus¬ 
ters from IWarren & Colei ( 2009 1. Although the authors origi¬ 
nally derived a linear calibration, the use of the updated [Fe/H] 
from iCarrera et al.l d2013l) for the open clusters is not compati¬ 
ble anymore with a single straight line. The difference between 
orange dots and red triangles, in Fig. |T0l looks more like a sys¬ 
tematic shift, rather than a change of slope. A similar shift can 
be observed in the [Ca/H] versus [Fe/H] plot shown in Fig. 11 
of ICarrera et alJ ( 2013 ). and it might be ascribed to inhomo¬ 
geneities in the [Fe/H] sources for open clusters, or to a sys¬ 
tematic effect due to age, compared to globular clusters. 

5. Conclusions 

Using RGB and RC stars from Baade’s window we have con¬ 
structed a CaT versus [Fe/H] calibration, valid in the metal rich 
regime up to [Fe/H ]=+0.7. Since bulge stars are rather sparse 
for [Fe/H]< -0.5, we used globular clusters to constrain the be¬ 
haviour of the calibration relation in the metal poor regime. We 
found that a single linear relation can be fitted across the whole 
range -2.3 <[Fe/H]< +0.7. A minor second order deviation 


from a single linear relation, however, gives a better fit, prop¬ 
erly taking into account the small difference between the linear 
best fit for globular clusters (equation Q} and for bulge giants 
(equation 0. It is worth emphasizing that both globular clus¬ 
ters and bulge stars have [Ca/Fe] ~ +0.3 in the common metal¬ 
licity range, therefore a possible dependence of the calibration 
upon the [Ca/Fe] ratio would not affect the combination of the 
two datasets. At higher metalliciti es the bulge [Ca/Fe] drops to 
zero (Gonzalez et alj2011 hL:2() 15 ) but the CaT-vs-Fe calibration 
found here is still almost linear, possibly due to the combined ef¬ 
fect of the strong saturation of CaT lines. The new calibration, 
represented by equation [3] allows one to derive the iron content 
of metal rich stars by means of low-resolution spectra. The latter 
are the only viable way to sample a large number of stars in high 
extinction regions, such as the Galactic bulge, or, in the future, 
in metal rich external spheroid such as the inner halo of M31. 
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Fig. 10. High resolution metallicity versus the reduced equiva¬ 
lent width for different star samples: bulge stars (black dots), 
globular (orange circles) and open (red triangles) clusters from 
[Warren & Cole ( 2009|), and globular clusters (blue stars) from 
Saviane et al. dSaviane et al.ll2012l) . Our CaT cali bration is plot¬ 
ted as a green line, while the cubic calibration from lSaviane et al] 
d2012l) is the dashed blue line. 
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